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Abstract 

Stability of neutrino masses and mixings with non-zero ^13 is an important issue in recent neu- 
trino physics. In context of this, a detail analysis on two- fold degenerate (m-i = 7712,7713) neutrino 

O 

^SJ ■ mass model is done with Renormaization Group Equation (RGE) in Minimal Supersymmetric 

^_i . 

C^ . Standard Model (MSSM). Considering the input values of neutrino masses and mixing angles at 

high scale predictions are studied at low energy scale using RGE. We have also found that two-fold 

C^ ' degenerate model with ^13 = can produce consistent results at low energy with non-zero ^13 from 

' ^ 'i the breaking of a generic A4 model. 
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I. INTRODUCTION 

Results of recent neutrino oscillation experiments have provided concrete of evidence 
which confirms the existence of the non-zero yet tiny neutrino masses [1]. Recenet neutrinio 
oscillation experiments like T2K [2], Double ChooZ [3], Daya-Bay [4] and RENO [5] have not 
only made the earlier predictions for neutrino parameters more precise, but also predicted 
non-zero value of the reactor mixing angle 6*13. The latest global fit value for 3a range of 
neutrino oscillation parameters [6] are as follows 

Aml^ = (7.00 - 8.09) x 10^^ 

Aml^{NH) = (2.27 - 2.69) x 10^^ 

AmlsilH) = (2.24 - 2.65) x lO'^ 

sin^^ia = 0.27 - 0.34 

sin2^23 = 0.34 - 0.67 

sin^^is = 0.016 - 0.030. (1) 

The above recent data have positive evidence for non-zero ^13 as well, which was earlier 
thought to be zero or negligibly small. The value of this mixing ^13 has serious impact in 
recent neutrino physics as it leads to the non-zero CP violating phase in the neutrino sector. 
It has also non-trivial impact on neutrino mass hierarchy as studied in recent paper [7]. 
Detailed analysis of this non-zero 613 using different approach is done viz. theoretical [8], 
phenomenological[9] before as well as after the discovery of this result in 2012. In view of 
the imortance of understanding the non-zero reactor mixing and CP violation in neutrino 
sector as mentioned above, present work is planned to present how a specific fi-r symmetric 
mass matrix [10] with two-fold degenerate mass eigenvalues can produce non-zero ^13 along 
with A?Ti2i through radiative corrections. This work is also emphasised that how a broken 
A4 symmetry at high scale can produce other neutrino oscillation parameters 6*23, ^12 and 
A?Ti23 at low energy scale through radiative corrections. 

This article is organized as follows: in section II we discuss briefly the broken A4 model 
which is used at high scale and about the RGE of mass eigen values and mixing to generate 
the non-zero reactor angle. In section III we discuss our numerical results and analysis. In 
the last section IV we outline the discussion and conclusions of our analysis. 



II. NEUTRINO MASS MODELS AND RGE 



A. A4 model for neutrino mass 



Type I seesaw framework is the simplest mechanism for generating tiny neutrino masses 
and mixing. The finite group of even permutation A4 can explain the fi — r symmetric mass 
matrix obtained from this type I seesaw mechanism. This group has 12 elements having 
4 irrudicible representations, with dimensions rij, such that '^^Ui = 12. The characters of 
4 representations are shown in table I. The complex number to is the cube root of unity. 
Under ^4 fermions and scalars of this model [11] transform as (z/j, /j) ~ 3, Zf ~ 1, V, 1", and 

TABLE I. Character table of A4 



Class x^'^ X^'^ X^') X^^^ 



Ci 1 1 1 3 



1 LO CO 







C2 

C3 1 a;2 w 

C4 1 1 1 -1 



3. Charged lepton mass matrix in this model can be written as [12]: 
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for Vi = V2 = V3 = V. Considering the charged lepton mass matrix for each allowed value, 
we have to consider the symmetry for M^. Let 
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f a-\- boj + cbj^ d 

e d a + bu"^ + cuj 



\ 



(3) 



be the Majorana neutrino mass matrix. Under A^ a, b, c, and {d,e,f) comes from 1, 1', 1" 
and 3 respectively. Now using the case b = c and e = / = [13] and by an ansatz a = 
and d = 3b we can arrived to the expression of Mj, as mentioned by [10] is 
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(4) 



which gives the two- fold degenerate mass eigenvalues {m,m,m'). 



B. RGE for neutrino mass eigenvalues and mixing 



The left-handed Majorana neutrino mass matrix miL which is generally obtained from 
seesaw mechanism at high scale Mr, is usually expressed in terms of K(t), the coefficient 
of the dimension five neutrino mass operator [14, 15] in a scale-dependent manner, t = 
\n{fi/lGeV) [16], 

mLLit)=vlKit), (5) 

where the vacuum expectation value (vev) is Vu = vq sin/3 and vq = 174 GeV in MSSM. The 
neutrino mass eigenvalues rrii and the PMNS mixing matrix Upmns [17] are then extracted 
through the diagonalisation of mLL{t) at every point in the energy scale t using equation 

(5), 

^diag ^ Diag(mi, 1712, m^) = VuLmLLVj^, (6) 

and the PMNS mixing matrix Upmns = Vl ^^ ^^^ basis where the charged lepton mass 
matrix is diagonal. The PMNS mixing matrix. 



Uel Ue2 Ue3 



u, 
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U^l U^2 U^3 

Url Ur2 Ur3 



(7) 



is usually parametrised in terms of the product of three rotations -R(^23), -^(^is) and R{6i2) 
(neglecting CP violating phases) by 
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U 



PMNS 



C13C12 C13S12 Sl3 

-C23S12 - C12S13S23 C12C23 - S12S13S23 C13S23 
■S12S23 ~ C12S13C23 — C12S23 — C23S13S12 C13C23 



\ 



(8) 



where Sij = sinOij and Qj = cos^ respectively. 

The RGEs for the eigenvalues of coefficient K(t) in equation (5), defined in the basis 
where the charged lepton mass matrix is diagonal, can be expressed as [18] 



d 



dt^^ = T^ E K-la! - Qd + QTrihl) + 2h}U%]K^. (9) 

Neglecting /i^ and h^ compared to h"^, and taking scale-independent vev as in equation (5), 
we have the complete RGEs for three neutrino mass eigenvalues, 

1^. = I^[(-^^? - 6^2 + 6/.?) + 2hlU'^,]m,. (10) 

The above equations together with evolution equations for mixing angles (19-20), are used 
for numerical analysis in the present work. 

The approximate analytical solution of equation (10) can be obtained by taking static 
mixing angle U^^ in the integration range, as [19] 

6 

niiito) = mi(tR)exp{-Igi + 6Ig2 - 6It)exp{-2U^Jr)- (H) 

The integrals in the above expression are usually defined as[16, 19] 

1 /•*« 

'to 

and. 



lAto) = Y^ I 9t{t)dt, (12) 



lAh) = ^^j\){t)dt, (13) 

where z = 1, 2, 3 and f = t,b,T respectively. For a two-fold degenerate neutrino masses , 
viz., rn^l^ = Diag.{m,ni,m') = Upmns''^llUpmns, equation (11) is further simplified to 
the following expressions 

mi (to) ^ m{tn){l + 25,(ci2Si3C23 - siaSss)^) + 0(5^), (14) 

m2(to) ~ -m{tn){l + 25,(c23Si3Si2 + C12S23)') + 0{S^), (15) 

m,ito) ^ m'(t^)(l + 25,(ci3C23)') + 0(5^). (16) 

While deriving the above expressions, the following approximations are made, 

exp(-2|f/^ip/^) ~ 1 - 2\Uri\'^Ir = 1 + 2\Uri\'^5r 



■'T1 



-5r = Ir^ —{m^/A-Kvf\n{MR/mt), 

cos^ p 

The quantity 5,- is always negative in MSSM and the approximation on 5t- taken here is 
vahd only if Iq is associated with the top-quark mass. From equations (14) and (15), the 
low energy solar neutrino mass scale is then obtained as 

Am^^(to) =ml-ml ^ 45^m^(cos26'i2(s^3 - 5^3033) + S13 sin 26^12 sin 26*23) + 0{5l). (17) 



C. Evolution equations for mixing angles 

The corresponding evolution equations for the PMNS matrix elements Ufi are given by 
|181 

at loTT^ ^-^, nik — mi 

where f = e, fi,T and i,k = 1,2,3 respectively. Here He is the Yukawa coupling matrices of 
the charged leptons in the diagonal basis and 

Neglecting h?^ and hi as before, and denoting A^i = ^^[^, , equation (18) simplifies to [18] 

Mci2[c2zSlzSi2UrlA^l - C23Sl3Cl2f/r2^32 + f^rlf^r2^2l], (19) 



dSv2_ _ 1 ,2, 

dt ~ 167r2 



ds 1 

-^ = T7r^^rC23C?3[Cl2f/rl^l + 512^7,2^32], (20) 

at lOTT^ 

^ = Y^^'4[-^12f/.l^l + Cl2f/.2A32]. (21) 

These equations are general and also independent of any neutrino mass model. 



D. Specific model for two- fold degenerate case 

We confine our analysis to a particular neutrino mass model [10] having two-fold degen- 
erate neutrino mass eigenvalues , viz., rrti = {m,m,m'). Then the master equation in (20) 
can be further simplified for input values m-i, = m' 7^ 0, S13 = 0, 

^ = Y^^rC23C?3[ci2Si2S23(A31 - A32)]. (22) 
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In both cases the mixing angle S13 evolves with energy scale, and this in turn induces solar 
mass scale through equation (17). 

III. NUMERICAL ANALYSIS AND RESULTS 

For a complete numerical analysis of the RGEs, equations (10),(19)-(21) for neutrino 
masses and mixing angles, we follow here two consecutive steps (i) bottom-up running [19] 
in the first place, and then (ii) top-down running [16] in the next. In the first step (i), 
the running of the RGEs for the third family Yukawa couplings {hf, hh, h^-) and three gauge 
couplings (91,92,93) in MSSM , are carried out from top-quark mass scale (to = Inm^) at 
low energy end to high energy scale M/j [19, 20]. In the present analysis we consider the high 
scale value as the unification scale M/j = 1.6 x lO^^GeV, with different tan/3 input values 
to check the stability of the model at low energy scale. For simplicity of the calculation the 
SUSY breaking scale is taken at the top-quark mass scale to = Inrnj [16, 19]. We adopt 
the standard procedure to get the values of gauge couplings at top-quark mass scale from 
the experimental CERN-LEP measurements at Mz, using one-loop RGEs, assuming the 
existence of a one-light Higgs doublet and five quark flavours below nit scale [19, 20]. Using 
CERN-LEP data, Mz = 91. 187 GeV, a s{Mz) = 0.118 ± 0.004, ai{Mz) = 127.9 ± 0.1, 
sin'^OwiMz) = 0.2316 ± 0.0003, and SM relations, 

1 3{l-sin^ewiMz)) 1 sin^OwiMz) , 2 



and,9i = Airai, (23) 



ai{Mz) 5 a{Mz) ' a2{Mz) a{Mz) 

we calculate the gauge couplings at Mz scale, ai(M^) = 0.0169586, 02 (M^) = 0.0337591, 
a^i^Mz) = 0.118. As already mentioned, we consider the existence of one light Higgs doublet 
[uh = 1) and five quark flavours {np = 5) in the scale Mz — rrit. Using one- loop RGEs of 
gauge couphngs, get 9i{mt) = 0.463751, 92{jnt) = 0.6513289 and (73(772^) = 1.1891996. Sim- 
ilarly, the Yukawa couplings are also evaluated at top-quark mass scale for input values of 
mt{mt) = 174 GeV, mb{mb) = 4.25 GeV, mt{mr) = 1.785GeV and the QED-QCD rescaling 
factors rjb = 1.55, rjr = 1.015 in the standard fashion [20], 



mt{mt)^Jl + tan^P 

^*^"^*) = — 174^;^^ — ■ 



hbimt) 



hr{mt) 



nibimt) ^1 + tarfi(3 
174 ' 

nirirnt) ^Jl + t(m?f3 



174 



(24) 



where mbirrit) 



rriiim,,) 



mAmt) 



mT(mT-) 



The one-loop renormahzation group equations 



% ' ' \ "/ rir 

(RGEs) for top quark, bottom quark and r-lepton Yukawa couphngs in the MSSM in the 



range of mass scales rrit < fj, < Mr are given by 
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(27) 



(28) 



The two-loop RGEs for the gauge couplings are similarly expressed in the range of mass 
scales rrit < 11 < Mr as 

d ■ 1 ^ ^ 



-T.9i 



dt" 167r2 
faa\ ( 



levr^ 



(29) 



i=i 



K 



6.6 

1 



\ 



/ 



j='t,b,T 
( 



and, aij 



(30) 



where 

/ \ / 

7.9 5.4 17 

,bij = 1.8 25 24 
2.2 9 14 

value of ht, hf,, h^, gi, g2, gz evaluated at high scale Mr = 1.6 x 10^^ from equation (25)-(27) 
and (29) are 



5.2 2.8 3.6 

6 6 2 

4 4 



\ 



/ 



ht{MR) = 0.977791667, /ife(MR) = 1.03416789, 
ht{MR) = 1.14931369, ^i(M^) = 0.718215168, 
gt{MR) = 0.726091146, ^3(Mi^) = 0.725136459. 
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In the second step (ii), the runnings of three neutrino masses {mi, 7712, ??23)and mixing angles 
(■512, ■523, ■513) are carried out together with the running of Yukawa and gauge couphngs, from 
high scale tij(= InMn) to low scale to- In this case, we use the input values of Yukawa and 
gauge couplings evaluated earlier at scale tji from the first stage running of RGEs in case 
(i). In principle , one can evaluate neutrino masses and mixing angles at every point of 
the energy scale. The same analysis is repeated for different values of tan/3 as well as the 
input value of high scale M/j. It can be noted that in the present problem, the running of 
other SUSY parameters such as Mq, M1/2, /U, are not required and hence, it is not necessary 
to supply their input values. As emphasised earlier, we study here a two-fold degenerate 
neutrino mass model , rrii = {■m,m,'m')aX high scale. 

We are now interested here to study radiative generation of both Am^i and sin 61^ for the 
case when m' ^ and S13 = 0. Such studies can give the possible origin of the solar scale 
as well as the reactor angle with broken A4 model. During the runnings of mass egenvalues 
and mixing angles in the RGEs in the present analysis, the non-zero input value of m' will 
first induce radiatively a non-zero values of S13 which in turn generates non-zero value of 
Am2i while running from high to low scale. 

We compare our numerical results with recent neutrino oscillation data mentioned in 
section I and observed the following results: 

• Variation of sin'^623 with energy shows that sin^623 is stable at low energy for different 
tan/3 values from 10 to 55 (figure 1). 

• Variation of sin^Ou with energy shows that sin'^612 is stable at low energy for different 
tan/3 values from 40 to 55 (figure 1). 

• Variation of Arn^^ with energy shows that Am'^^ is stable at low energy for different 
tan/3 values from 10 to 52 (figure 1). 

• Generation of Am2i with energy shows that Arn^i is stable at low energy for the tan/3 
values 50 to 55 (figure 1). 

• Generation of sin^di^ with energy shows that sin^Oi^ is stable at low energy for the 
tan/3= 55 (figure 2). 

It is observed from the figure that all the parametrs stable under radiative corrections with 
taUyS = 55 except the fact that prediction for Arri^^ are slightly deviated from 3cr range of 




FIG. 1. Radiative generation of the sin?6i2, si'n?92z, ■^'n-li 'ind ATTi^a at low energy for different 
values of tan/3 




0.372 0.374 0.376 0.3,78 0.38 0.382 0.384 



FIG. 2. Radiative generation of sin^Oi^ for different tan/3 values at low energy, and variation of 
siv?6i^ with respect to sin?'6i2, sin^023 and /S.m^i for tan/3=55 

experimental data. Therefore, it is very likely that more precise future data from neutrino 
oscillation experiments can expain these exceptions in the present analysis. 
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IV. DISCUSSION 

Analysis of the effect of neutrino mass eigenvalues and mixings in case of two-fold de- 
generate neutrinos are done with RGE in the broken A^ model. It is observed that all the 
parameters are stable at low slcale under radiative corrections for tan/3 around 55 except 
the fact that for Amga the prediction are slightly below the 3a range of experimental data. 
Apart from this exception, all other predicted values of the neutrino parameters are consis- 
tent with neutrino oscillation data. In view of above, the scenario of two-fold degenerate 
neutrinos can survive in nature within the framework of Type I seesaw mechanism. However, 
here we have not made any attempt to explain the CP violating phase. As an extension of 
this work, one can incorporate the RGE for phases which we have left for future studies. 
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